Mesoscopic fluctuations in superconducting dots at finite temperatures 
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We study the thermodynamics of ultrasmaU metaUic grains with the mean level spacing 5 compa- 
rable or larger than the pairing correlation energy in the whole range of temperatures. A complete 
picture of the thermodynamics in such systems is given taking into account the effects of disorder, 
parity and classical and quantum fluctuations. Both spin susceptibility and specific heat turn out 
to be sensitive probes to detect superconducting correlations in such samples. 



Superconducting coherence is expected to be com- 
pletely washed out in small metallic grains, where the 
average level spacing S ^ 1/{JV{0)V) may be compa- 
rable or even larger than the BCS energy scale A [Q. 
Despite of that, a series of experiments by Ralph, 
Black and Tinkham (RET) revealed a gap related to su- 
perconductivity in individual nanosized metallic grains, 
(5 < A, with large electrostatic energy fixing the total 
number N of electrons. The BCS description of super- 
conducting coherence is no longer valid and the charac- 
terization of "superconductivity" from the bulk to the 
ultrasmall-grain regime has been studied in recent theo- 
retical works J0-01 . Starting from the standard pairing 
Hamiltonian |M, several spectral features at T = were 
calculated, which turn out to be universal functions of 
the single scaling parameter (5/A ft is remarkable 
that pairing determines strong fluctuational "supercon- 
ductivity" even in ultrasmaU grains Q, where S > A 
(we call this the "dot" regime). As in micron-size sam- 
ples parity effects are present, i.e. physical properties 
depend on N being odd or even. An exact solution for 
the pairing problem in finite systems was found long ago 
by Richardson and Sherman (RS) ||]. 

In this work the thermodynamics of an ensemble of 
monodispersed superconducting dots is studied. The spe- 
cific heat cv{T) (see Fig. P and the spin susceptibility 
x{T) (see Fig. ^ are obtained in the whole tempera- 
ture range. Our work is motivated by the fact that ther- 
modynamic properties are a unique experimental tool in 
detecting unambiguous traces of superconducting corre- 
lations in the dot regime, where tunneling spectroscopy 
on individual dots is not sensitive enough . Di Lorenzo 
et al. [0 proposed that a possible signature of pairing 
for (5 > A, is the reentrant behavior of the spin suscep- 
tibility Xo{T) in dots with odd N, which is due to the 
combined parity and interaction effects, using a simple 
equally spaced spectrum for the single particle energies. 

In ensembles of grains the single-particle spectrum is 
statistically distributed and characterized by level repul- 
sion ||T^. The statistics has a universal character due to 
general symmetry properties and is described by Random 
Matrix Theory (RMT). The thermodynamics of ensem- 



bles of normal metal grains was studied long ago 
but only in recent years the first experimental evidence 
of mesoscopic fluctuations and level repulsion in normal 
metal grains was provided with the observation of the 
T^-dependence of the specific heat and of other thermo- 
dynamic quantities Pioneering experimental works 
on thermal properties of small superconducting grains 
date back to the 80s . 

Concerning superconductors, level statistics was found 
to yield enhanced superconductivity and parity effects 
in micron-size samples [^5[. Similar results have been 
recently found also in the dot regime |]l6| at T = 0. Here 
we focus on ensembles of dots at finite temperature. They 
exhibit detectable features characteristic of the interplay 
between pairing interaction and the universal statistics of 
mesoscopic fluctuations, which determines the physics. 

Our starting point is the Hamiltonian 

7i = ^(cq- cr/i_Bi/)c]j_^Ca,^- A rt T , (1) 

a, a 

where a spans a shell of pairs (cr — ±) of single par- 
ticle energy levels with statistically distributed energies 

and annihilation operator Ca.a- The magnetic field 
H enters in the Zeeman form (pB is the Bohr magne- 
ton). For H = the pair of states a — ± are degenerate 
and time-reversed. The operator T — J2a=i Ca,-Ca,+ ap- 
pears in the interaction term, which scatters pairs with 
amplitude A, so a pair of levels cr = ± occupied by 
a single electron is blocked Relevant energy scales 
are the average level spacing 5 and the BCS energy, de- 
fined as A = ri(5/[2 sinh((5/A)] . This model stems from 
a rather general description of electron-electron interac- 
tion in the metallic regime (large dimensionless conduc- 
tance, g — Et/S :s> 1, where Et is the Thouless energy) 
which has been recently proposed . The Hamiltonian 
Eq.(^ describes universal properties of pairing interac- 
tion in disordered metallic dots, with fixed N and negli- 
gible exchange interaction. We assume that the energies 

are distributed according to the Gaussian Orthogo- 
nal Ensemble (GOE), which describes mesoscopic fluc- 
tuations when no spin-orbit interaction is present and 
time-reversal symmetry is preserved pd| . 
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Level statistics has a drastic effect if T <C 5, where the 
thermodynamics is determined by samples with Si, the 
closest level spacing to the Fermi energy, much smaller 
than the average S. The relevant excitations involve 
only few levels around the Fermi energy, and excitation 
energies can be estimated by solving the problem in a 
small shell, 51' ^ fi, with renormalized coupling con- 
stant X = Xn ^ Ao' P,|6|,p^ . Let's consider dots with 
even N and estimate the leading low temperature con- 
tribution to cve{T). In the regime (5 ^ A, the system 
is first renormalized to an effective one with two elec- 
trons in two doubly degenerate levels, and interaction 
A2. We have to consider six states, and the relative ex- 
citation energies: ground state, excited state of paired 
electrons {Ep{Si) = 2{Si + ^2)^^'^) and four states corre- 
sponding to a broken pair with different spin configura- 
tions (EhpiSi) = A2 -I- Ep{6i)/2). Notice that even if the 
coupling is weak, A2 <IC 5, dominant configurations have 
Si ^ X2 ■ The response of a set of grains is evaluated by 
retaining only fluctuations of Si, governed by the level 
spacing distribution Pj^ (Si). In this way we obtain the 
analytic form for T <C (5. The leading term for T ^ A2 
reads 

cye(T) = Ae-#^ (2) 

and shows a gap 2A2. In the limit X2 ^ T <^ S the full 
analytic result reproduces the known result for normal 
metal grains For cvo{T) we have to consider an 

effective three (fluctuating) level problem, and we obtain 

cvoiT) = ^.2c(3)(^)' + o(^)', (3) 

which shows no trace of a gap. The contribution of su- 
perconducting correlations for T ^ A3 is a correction of 
the result for normal grains ||l^ , because the most im- 
portant low-energy excitations are obtained by moving 
the unpaired electron. 

Dots at intermediate temperatures can be studied nu- 
merically by the massive use of the RS exact solution. 
Noninteracting random spectra are obtained by diagonal- 
izing 500 x 500 random matrices belonging to the GOE 
and taking the central levels of the resulting semicircular 
distribution, in order to prevent undesired border effects. 
For each disorder realization we evaluate the partition 
function from the universal Hamiltonian Eq. (0) using the 
RS solution and calculate the average free energy. For the 
curve S/A — 50 shown in Fig. ^we used sets of 50 levels. 
At intermediate temperatures (T < 0.25(5 for the curve 
in Fig. ^ we used 500 realizations and an energy cutoff 
5S. At larger coupling A, corresponding to the crossover 
region 6/ A > 1, and when excited states with random 
spectra are concerned, singularities in the RS equations 
become untractable. Then the low-temperature behav- 
ior for S/A > 1 shown in Fig. Il| was studied following 



Ref. for each realization of disorder and configura- 
tion of blocked levels the interaction is scaled down to 
some low-energy cutoff and the effective small (< 12 lev- 
els) system, and is then diagonalized. 

Both methods described above cannot be used to ob- 
tain results at temperatures T ~ 5, A where cv (T) is ex- 
pected to show anomalies. In this regime a huge number 
of states (at least 5 • lO'' already at T ~ ^) are needed for 
each single disorder realization. Then for high tempera- 
ture we use a functional technique in the parity-projected 
(PP) grand canonical ensemble. Following Janko et al. |^ 
we introduce the even and odd N partition function 

00 

Ze/o{T, m) - 2 E ^'""^^^^ ± e"^]Z(T, N) , (4) 

for each given realization of disorder. Here Z{T, N) is the 
canonical partition function so Ze/o{T,^) involve only 
sectors where N is not fixed but its parity is. The effec- 
tive imaginary-time action is obtained from Eq. (^ in the 
standard way and Ze/o(T, /i) can be expressed as a path 
integrals over a Hubbard-Stratonovic (HS) auxiliary field 
A(t) [0. This non- linear functional integral is evaluated 
in the RPA', an approximation which has been widely 
employed to study the Anderson model |lj]. We extract 
the static part of the HS field in the Matsubara-Fourier 
representation, A(t) = Aq + ^„_^q A(co'„) exp(— ia;„T) , 
evaluate quantum corrections by including the contribu- 
tion of small-amplitude deviations around a generic static 
path Aq and perform numerically the remaining ordi- 
nary integral over Aq. We used systems of 100 levels 
distributed according to GOE. Finally the resulting free 
energy is averaged over 400 samples. Results for cve{T) 
are shown in Fig.^. For large enough - albeit still nano- 
sized - grains a finite temperature anomalous enhance- 
ment develops, which can be attributed to the presence 
of pairing interactions. 

We now study the spin susceptibility. Using the same 
arguments leading to Eqs.(^,||) we can conclude that at 
low temperatures the susceptibility of grains with even N 
is exponentially suppressed, so we concentrate on Xo{T) ■ 
For equally spaced single particle spectrum Xo{T) shows 
a reentrant behavior for non-vanishing pairing interac- 
tion, because at low temperature the leading contribu- 
tion is the 1/T Curie-like coming from the unpaired elec- 
tron, whereas the response of the paired electrons ex- 
ponentially increases with T and becomes dominant at 
larger temperatures |l^. As compared with results for 
equally spaced spectra, we notice that level statistics trig- 
gers two competing effects: the availability of low-lying 
excitations tends to increase Xo{T) , which would wash 
out the reentrance; on the other hand disorder enforces 
superconductivity, which would determine the opposite 
trend. 

We evaluate Xo{T) using the PP partition function 
Eq.(|) in the RPA'. Results (Fig. |[) clearly show that 
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Xo{T) is reentrant also in disordered samples. The reen- 
trance is slightly less pronounced than for regular spec- 
tra and it is present always but for unphysically small 
grains. Notice that the "superconducting" Xo{T) is well 
suppressed with respect to normal metal grains, confirm- 
ing the expectation that disorder favors superconductiv- 
ity. The leading behavior of Xo{T) for T <C (5 is unaf- 
fected by level statistics. 

We now discuss experimental signatures of supercon- 
ducting correlations. The main point is that in ensembles 
of dots available in experiments roughly 50% have odd 
iV and 50% have even N. At low temperatures the effect 
of pairing is detectable in the total cyiT): in an ensem- 
ble of normal metal grains the relative contributions to 
the ensemble specific-heat of dots with even N is about 
63% Pairing suppresses this contribution, Eq.(2), 

whereas in practice it does not affect cvo{T), Eq.(3). 
Thus pairing reduces the total cy(T) of ^ 2/3. An- 
other signature is that cv(T) increases with an applied 
magnetic field. Other signatures can be found at inter- 
mediate T: The total susceptibility xC^) is expected to 
show reentrant behavior (x(T) « Xo{T), since XeiT) is 
exponentially suppressed), except for extremely weak in- 
teractions or extremely small grains. Another signature 
is the increase of x{T) with the applied magnetic field. 
In the same temperature regime ensembles of dots with 
show an anomalous enhancement arising from 
the contribution of cve{T) which can be shifted towards 
lower temperatures and the progressively suppressed by 
applying a magnetic field. The whole phenomenology 
could provide unambiguous evidence of the presence of 
pairing correlations even in the dot regime. 

The procedure used to obtain the analytic results for 
T <C (5, neglects fluctuations of the renormalized cou- 
pling A2 (A3 for dots with odd N). They arise from 
fluctuations of larger shells of levels which may be sys - 
tematically accounted for. In normal metal grains [[l2[ , 
they can be neglected for T <^ 5. We checked numer- 
ically that cve is exponentially suppressed and that for 
(5/A ^ I the exponent roughly agrees with the scaling 
form discussed in Ref . . We used a large number of re- 
alizations lO"' — 10^) since at very low T, the system 
is very sensitive to the statistic. Detailed account will be 
given elsewhere. 

Some remarks on the results presented above are 
clearly stated if we draw a comparison with a dot having 
equally spaced noninteracting levels, ta = oi5. In this 
case all the thermodynamic quantities show gapped be- 
havior for T <^ 5. This is essentially due to the finite 
level spacing Jf^,|l^, for instance the even specific heat 
has a gap 5 + X2 ■ On the contrary, in ensembles of grains 
gapped behavior is not the rule. When a gap is present, 
it is due only to the pairing interaction, as in the case of 

Cye(T), Eq.(|). 

In Fig.(^) we present results at finite temperatures for 
cve{T) and Xo{T) for regular spectra, = a5. In this 



case the exact RS solution allows to reach temperatures 
slightly larger than 5] larger temperatures (T > 5) were 
studied by the PP approach, Eq. In the inset we 

show Xo{T), calculated using the canonical RS solution, 
and the PP grand canonical approach in the RPA' and 
in the Static Path Approximation (SPA) [|To|. In the 
SPA fluctuations around the static path Aq are neglected 
when evaluating the PP partition function Eq. (jij). The 
SPA accounts for thermal fluctuations and it has been 
successfully employed to study the thermodynamics of 
micrometer size grains p9| . By comparing the SPA with 
the RPA' results we can quantify quantum corrections. 
As expected quantum superconducting fluctuations tend 
to enhance the features of pairing: the reentrance in 
Xo{T) is deepened compared to the SPA result getting 
closer to the exact result for T ^ 5. Differences be- 
tween the RPA' and the exact result reflect the different 
sets of excitations available in the canonical and in the 
PP grand-canonical ensembles. The discrepancy is small 
and probably undetectable when pairing interaction is 
active, so RS solution and RPA' can be used in combina- 
tion. Small differences are present also in cve (T) (Fig.^) 
and they are expected to persist at T ^ ^. In this latter 
regime cve{T) approaches the behavior of normal met- 
als, where it is known that the canonical cv{T) is smaller 
than the grand-canonical one, asymptotically by a quan- 
tity ifcB @. 

We now focus on the anomaly at intermediate temper- 
atures appearing in cve{T). The qualitative behavior of 
cve{T) for (5/ A ^ 1 is reminiscent of what observed in 
bulk superconductors, but this is partly due to the special 
choice of the equally spaced noninteracting spectrum. In 
fact the T < 6 enhancement in the specific heat is present 
also if the pairing interaction is switched off. This artifact 
is not present for normal grains with levels Eq, distributed 
according to the GOE statistics (symbols in FigJ|) [p^ . 
Thus the analysis of the effects of level statistics is neces- 
sary for a reliable study of the thermodynamics of small 
grains and gives an unequivocal discrimination of the ef- 
fects of the pairing interaction. 

In summary we studied the thermodynamics of ensem- 
bles of ultrasmall superconducting grains, addressing for 
the first time this problem in the dot regime. We draw 
a picture of the physics which comes from the interplay 
of disorder, strong pairing fiuctuations and fixed num- 
ber of electrons in each dot. Thermodynamic proper- 
ties are a unique experimental tool in detecting unam- 
biguous traces of superconducting correlations in dots. 
We found that several signatures of pairing (the low and 
intermediate temperature anomalies of the specific heat 
and the reentrant spin susceptibility) are experimentally 
detectable in ensembles of dots, despite superconduc- 
tivity in the BCS sense breaks down. Finally the low- 
temperature behavior of cy(T) may provide direct evi- 
dence of the effects of level statistics in interacting 
metallic systems. 
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FIG. 1. Rescaled specific heat cve5/'^T vs. temperature in 
units of 5, disordered case. Circles: the RPA' results; trian- 
gles: the exact canonical results. 
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FIG. 2. Spin susceptibility for dots with odd N . Equally 
spaced noninteracting single-particle levels case as given by 
the RS canonical solution (dashed lines), and the RPA' 
(solid), axe compared with results for GOE distributed sin- 
gle-particle levels obtained by RPA' (circles) and with the 
exact canonical calculation (triangles). 




FIG. 3. Thermodynamic quantities for dots with equally 
spaced spectra, calculated by the RS solution (dashed lines), 
the SPA (dotted) and the RPA' (solid) as a function of tem- 
perature in units of 5. Main figure: the rescaled specific heat 
cve/'yT. The curve for normal metal disordered dots is plot- 
ted for comparison (triangles), see text. In the inset the odd 
spin susceptibility Xo{T), see text. 
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